Nonstoichiometric doping and Bi antisite defect in single crystal Bi 2 Se 3 
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We studied the defects of Bi2Se3 generated from Bridgman growth of stoichiometric and nonsto- 
ichiometric self-fluxes. Growth habit, lattice size, and transport properties are strongly affected by 
the types of defect generated. Major defect types of Bis e antisite and partial Bi2-layer intercalation 
are identified through combined studies of direct atomic-scale imaging with scanning transmission 
electron microscopy (STEM) in conjunction with energy-dispersive X-ray spectroscopy (STEM- 
EDX), X-ray diffraction, and Hall effect measurements. We propose a consistent explanation to the 
origin of defect type, growth morphology, and transport property. 



PACS numbers: 61.72.Ff, 74.25.F-, 74.70.Xa 

Bi 2 Se3 is the second-generation topological insulator 
with a nearly idealized single Dirac cone.pQ Bi2Se3 is 
also the sister compound of Bi2Te3 that has been demon- 
strated to show high thermoelectric figure-of-merit (ZT) 
for commercial applications. [2] Both Bi2Se3 and Bi2Te3 
samples show nearly identical crystal structure composed 
of Bi-(Se,Te) quintuple layers with a van der Waals gap 
between Se(Te)l, as shown in Fig. []Ja). While Bi 2 Te 3 
has been shown with n- or p-type conduction through 
fine-tuning of nonstoichiometry of Bi/Te ratio, i.e., to 
dope the system by generating antisite defects of Bi^e 
(p-type) or Te^ (n- type), [3 all Bi 2 Se3 samples reported 
since 1960 showed n-type semiconducting behavior with- 
out exception, and the lowest n-type carrier concentra- 
tion reported is at the order of 10 16 -10 17 cm -3 . [4 

Topological insulator is expected to be a band insu- 
lator with surface conduction and spin-polarized surface 
state. pQ When the Fermi level is higher in the conduc- 
tion band as a result of intrinsic n-type doping, it be- 
comes difficult on data interpretation for results contain- 
ing contributions from both the bulk and the surface. An- 
gle resolved photoemission spectroscopy (ARPES) stud- 
ies require n-type doping to map out the shape of the 
Dirac cone for Bi2Se3, but transport property measure- 
ment faces more challenge on distinguishing conductivi- 
ties from the bulk and the surface. Additional hole dop- 
ing methods through Ca substitution and NO2 oxidation 
have been applied to provide the necessary fine-tuning of 
Fermi level near the Dirac cone. [5j [6] It is desirable to 
have complete carrier type and level control for the topo- 
logical insulator Bi2Se3 at the stage of crystal growth. 

In an effort to grow Bi2Se3 single crystals of controlled 
carrier type and level for surface- sensitive studies, and 
to clarify the nature of defects responsible for the per- 
sistent n-type character, we have studied single crystals 
grown using the Bridgman technique with intentionally 



generated defects through non-stoichiometric self- fluxes. 
With the integrated studies of direct defect observation 
through STEM-EDX imaging, X-ray diffraction lattice 
size analysis, and transport property investigation, we 
were able to find reasonable correlations between the ob- 
served physical properties and the distribution of defects 
in the form of either Bis e antisite or partially intercalated 
Bi2-layer patches. 

Single crystals studied in this work were grown us- 
ing the vertical Bridgman method. The initial com- 
pounds were high purity bismuth (99.999%) and sele- 
nium (99.999%) powder mixed thoroughly in an argon- 
filled glove box. Three groups of sample were prepared 
in different Bi:Se molar ratios, including Bi:Se molar ra- 
tios of Bi:Se=2:3, (2+x):3, and 2:(3+x) with x=0.05-0.2. 
The mixed powder were sealed in quartz tubes under 
vacuum of about 3xl0 -2 torr after multiple argon gas 
purging cycles, pre-reacted at 650° C for 18 hours in a 
box furnace, and furnace-cooled to room temperature. 
Single crystals were grown with a vertical Bridgman fur- 
nace starting from the pre-reacted powder and vacuum- 
sealed in quartz tubes of 10 cm long and 1.6 cm inner 
diameter. The temperature profile of the Bridgman fur- 
nace used for the whole series was maintained at 850° C- 
650° C in 25 cm. Initial complete melting was achieved 
at 850° C for 24 hours to ensure complete reaction and 
mixing. The temperature gradient of l°C/cm was pro- 
grammed around the solidification point near 705° C, and 
the quartz tube was then slowly lowered into the cooling 
zone at a rate of ~0.5 mm/h. Crystal structure was ex- 
amined by X-ray diffraction with in-house powder X-ray 
diffractometer (Bruker D8). Transport properties were 
examined using a standard four-probe lock-in technique 
with magnetic field up to 2 Tesla and temperature down 
to 2 K. High-angle annular dark- field (HAADF) imaging 
and chemical mapping at atomic-column resolution were 
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FIG. 1. (Color online) (a) Crystal structure of Bi 2 Se3 can be 
described with space group R3m, where quintuple layer is the 
building block with Se2(blue) in the middle and Sel (green) 
near the van der Waals gap. The inset shows the BiSe6 oc- 
tahedron and the Sel (or Bisei antisite in dashed circle) en- 
vironment, (b) Crystals grown using Bi:Se=2:3 and 2:(3+x) 
of stoichiometric and Se-rich fluxes, with (001) plane index 
mostly perpendicular to the growth direction, (c) Crystals 
grown using Bi:Se=(2+x):3 of Bi-rich flux show (001) plane 
index along the growth direction. 



performed, with JEOL-2100F microscope equipped with 
a probe Cs corrector in conjunction with EDX as detailed 
previously. |5] . Chemical analysis was performed using 
Electron probe microanalysis (EPMA) on pieces sepa- 
rated from the as-grown cone-shaped solidified boule. It 
is interesting to note that the Se-rich flux growth shows 
growth morphology mostly, but not always, with (001) 
plane index perpendicular to the growth direction, as 
shown in Fig. Qb). On the other hand, the growth from 
Bi-rich flux has (001) plane index parallel to the growth 
direction always and can be cleaved easily into parallel 
discs, as shown in Fig. [3c). 

STEM-HAADF images for samples selected from dif- 
ferent initial growth flux ratios and positions are summa- 
rized in Fig. [2] Fig. ^a) can be viewed as the [100] pro- 
jection sequenced as quintuple Sel-Bi-Se2-Bi-Sel, which 
shows two bright spots corresponding to the heavier Bi 
atoms and the three weaker ones corresponding to the 
Se atoms as expected, as well as the van der Waals gaps 
in between the quintuple layers. For crystals grown from 
Bi-rich flux in lower magnification, as shown in Fig.[2|b), 
bright domain boundaries can be found and verified to 
be Bi, using STEM-EDX characterizations. As expected 
for the Bi-rich flux growth, the excess Bi may aggregate 
into grain boundaries during cooling, most probably due 
to the low melting point of Bi (271° C), while crystal so- 
lidification occurs near 705° C. Moreover, low density of 




FIG. 2. (Color online) STEM-HAADF images of Bi 2 Se 3 . (a) 
High resolution image shows that the quintuple of Bi2Se3 con- 
sists of two bright spots (Bi) and three fainter spots (Se). (b) 
For the Bi-rich flux growth, high density of Bi can be found 
in the grain boundary as shown by the thick and bright lines 
in low magnification, (c) Low density of Bi2-layer interca- 
lated patches can always be found in the van der Waals gap 
for the as-grown crystals, whether using Bi- or Se-rich flux. 
The intercalated Bi2-layer patches would deform the quintu- 
ple layer locally by opening up the van der Waals gap. (d) 
Large area of perfectly ordered quintuple layers without inter- 
calated Bi2-layer patches induced deformation can be found 
in samples after 300° C annealing. 



Bi2-layer patches were found in the van der Waals gap 
sporadically, for both the Bi-rich and Se-rich flux grown 
crystals, as shown in Fig. [2^c) . To accommodate the in- 
tercalated Bi2-layer patches, the quintuple units near the 
enlarged van der Waals gap are locally distorted. Inter- 
estingly, large areas of structure free from imperfections 
could only be found after 300° C annealing for both the 
Se-rich and Bi-rich grown crystals, as shown in Fig.[2|d), 
which indicates that both Bi excess in the grain bound- 
aries and the intercalated Bi2-layer patches had been re- 
moved accordingly. 

Lattice parameters of c-axis have been examined us- 
ing X-ray diffraction for pieces selected from different 
sections of the grown boule. Additional annealing pro- 
cess has been applied at 300° C to 600° C in an evacuated 
sealed tube for part of the as-grown crystals. Lattice 
parameters are summarized in Fig. [3] and plotted in the 
a- and c-axes against lattice volume. Interestingly, all 
studied samples are clustered in groups of Bi-rich flux 
growth before annealing (Bi-rich), Se-rich flux growth be- 
fore annealing (Se-rich), and Se-rich plus 300-600° C post 
annealing (Se-rich+300-600°C), in a trend of increasing 
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FIG. 3. (Color online) Lattice parameters vs. lattice volume 
of all crystals studied, where c-axes are in empty circles and a- 
axes are in empty squares. All studied samples are clustered 
in groups of Bi-rich flux growth before annealing (Bi-rich), 
Se-rich flux growth before annealing (Se-rich), and Se-rich 
plus post annealing at 300-600° C (Se-rich+300-600°C), in an 
increasing trend as indicated by the connected experimental 
data points. A list of Bi2Se3 lattice parameters reported in 
the literature since 1960 are shown in various solid symbols. [9] 



volume and c-axis. In addition, a collection of lattice pa- 
rameters of Bi 2 Se3 published since 1960 are displayed to- 
gether and found falling nicely along the same increasing 
trend. [9] It is clear that many of the inconsistent results 
published previously could be due to the subtle differ- 
ences in the nonstoichiometry of Bi and Se. 

We may postulate the reason why crystals from Bi- 
rich flux growth have longer c-axes is due to the larger 
amount of Bi intercalated into the van der Waals gaps, 
mostly in the form of randomly distributed patches of 
neutral metal Bi 2 -layer as revealed in the STEM pic- 
ture of Fig. [5|c), but not enough to form the metastable 
phases of staged (Bi2) m (E>i2Se3) n at the high solidifica- 
tion point near 705° C. [10] It should be noted that the 
staged phases beyond Bi2Se3 can only be prepared at the 
much lower temperature range of ^-400-560° C as interca- 
lated metastable compounds. [TT] In particular, the aver- 
age lattice volume increase should be viewed as randomly 
distributed intercalated Bi2-layer patches without chang- 
ing the original R3m symmetry, especially when the local 
dilation of van der Waals gap occurs only near regions 
with intercalated Bi2-layer patches. In addition, 300° C 
annealing reduces the lattice size significantly, which in- 
dicates that a small amount of intercalated Bi can be 
removed with a temperature above the melting point of 
Bi effectively, as supported also by the creation of large 
area of defect-free structure after 300° C annealing shown 
in Fig. [2^d). On the other hand, high annealing tem- 
perature up to 600° C starts to decompose the sample 
partially, presumably as a result of too many Se vacan- 
cies created in an evacuated environment, which has also 
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FIG. 4. (Color online) Drude mobility (p) vs. carrier concen- 
tration (n) for all samples studied are shown; all are n-type. 
Crystals grown from Bi-rich flux (Bi-rich) have the highest 
carrier concentration but the lowest mobility, and post 300° C 
annealing (Bi-rich+300°C) reduces the carrier concentrations 
to the same level of ^-15xl0 18 . Crystals from stoichiometric 
and Se-rich flux growth (Se-rich) have the lowest carrier con- 
centration and highest mobility, and the 300° C post annealing 
(Se3. io+300 C) shows no significant impact on modifying the 
carrier concentration and mobility. 



been verified by thermogravimetric analysis under argon 
flow with the following X-ray structure analysis. 

All single crystal samples show metallic behavior down 
to about 20 K (not shown, similar to those published 
earlier) [12] , below which the resistivity p becomes nearly 
temperature-independent, as expected in the impurity 
scattering regime. Fig. [4] summarizes the relationship 
between carrier concentration n and mobility p for all 
samples studied. The carrier concentration n was deter- 
mined from Hall effect measurement at 2 K, and the cor- 
responding Drude mobility p was estimated using Drude 
formula 1/p = nep. For the stoichiometric and Se-rich 
flux growth, samples are clustered at the low carrier con- 
centration ~ -5 x 10 18 /cm 3 and high mobility regime. 
On the other hand, crystals grown from Bi-rich flux have 
the highest carrier concentrations and the lowest mobil- 
ity. In addition, 300° C annealing has no significant im- 
pact on the already low carrier concentration from Se-rich 
flux growth, while the carrier concentrations for Bi-rich 
flux growth crystals are reduced to a common lower level 
near ~ —15 x 10 18 /cm 3 . 

Based on evidences from partial decomposition after 
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600° C annealing in an evacuated tube, thermogravimet- 
ric analysis, and lattice size analysis, Se vacancy must 
be limited to a low threshold before the occurrence of 
decomposition. This is because Se has a vapor point 
of 685° C, below the solidification point of Bi2Se3 near 
705° C, which leads to the occurrence of solidification at 
the liquid-gas phase boundary in the Bridgman growth. 
It is nearly impossible to maintain a perfect Bi-Se stoi- 
chiometry upon Bridgman solidification near 705° C, es- 
pecially for the relatively weakly bonded Sel at the van 
der Waals gap. We conclude that the Se vacancy is im- 
possible to avoid under the current Bridgman growth 
condition of limited Se vapor pressure. The Se vacancy 
level must be determined mainly by the Se partial pres- 
sure during Bridgman growth, instead of the initial Bi 
to Se ratio, not to mention the fact that the liquid Bi in 
the flux prefers to intercalate into the van der Waals gap 
after the Bi2Se3 solidification, instead of forming Bi-Se 
quintuple layer of high Se vacancy level. Although lack- 
ing rigorous quantitative evidence, current results seem 
to agree with the commonly accepted assumption that 
n-type carriers are generated from Se- vacancy formation, 
i.e., each Se vacancy donates at most two itinerant elec- 
trons as Se^e ^V^ e +Se(3)+2e - in the Kroger- Vink no- 
tation of various doping efficiency. [13 

In addition to the finding of intercalated Bi2-layer 
patches, as revealed in Fig. |2|c), we also observed Bis e 
antisite defects in the bulk of ~25 nm thickness along 
the incident beam direction, as shown in Fig. [5j us- 
ing atomically-resolved STEM-EDX mapping. We no- 
tice that some columns of supposed Sel arrays near the 
van der Waals gap show unexpected enhanced contrast, 
as shown in Fig. |5|a). The atomically resolved chemi- 
cal mapping by STEM-EDX shown in the right panels 
of Fig. [5|a) reveals that the enhanced contrast at the 
Sel site actually corresponds to Bi dominant character, 
which provides a direct proof of partial substitution of Bi 
to the Sel site, i.e., the Bis- e antisite defects formation. 
The observed Sel columns with an enhanced contrast to- 
wards the sample edge (Fig.^b)) signify the existence of 
more antisite defects near the surface layers. The antisite 
defects both in the bulk and near the surface region has 
been a general phenomenon throughout our STEM-EDX 
characterization of various batches of samples. 

The importance of antisite defect in the tetradymite- 
type compounds has been discussed in the scope of nar- 
row band gap semiconductor previously, [13] yet quanti- 
tative evidence is difficult to obtain at the point defect 
level. Considering the number of outer shell electrons for 
Bi(6s 2 6p 3 ) and Se(4s 2 4p 4 ), simple argument borrowed 
from extrinsic semiconductor suggests that Bis e antisite 
can be an acceptor. However, it would be a picture too 
simplified when Bis e antisite defect is actually formed 
by the Bi substitution only to the Sel vacancy site near 
the van der Waals gap, i.e., Bis e i antisite, and chem- 
ically, the electronegativity difference between Se(2.55) 
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FIG. 5. (Color online) (a) STEM-HAADF image of a bulk 
region with thickness ^25 nm along the incident beam di- 
rection. Sel column pointed out by the white arrow shows 
an anomalously enhanced imaging contrast. Further STEM- 
EDX chemical mappings of the corresponding boxed region 
are enlarged in the right panels (cyan, Bi map; yellow, Se 
map). The gray (white) circles in the mapping represent the 
Bi (Se) columns. The enhanced contrasts of the Sel columns 
(signified by the red Sel text) are shown to be correlated with 
the presence of Bi on the sites, i.e., the formation of Bi,s e an- 
tisite defects, (b) The enhanced Sel contrasts are more regu- 
larly observed at the edge of the same specimen as indicated 
by the white arrows, with corresponding enhanced blow-up of 
the white boxed region in the inset in lower left corner. 

and Bi(2.02) is too large for Bi to accept electrons fur- 
ther from the neighboring Bi/Se atoms. Originally, Sel 
uses its four outer shell electrons in 4p 4 to bond with the 
neighboring Bi's and the two 4s 2 electrons form lone pairs 
exposed to the van der Waals gap, as shown in the inset 
of Fig. [TJa).[l4] The bonding nature of the substituted 
Bi in Bis e i antisite becomes unknown in view of a vir- 
tual octahedral environment of Bi(Bi3Sels) as depicted 
in the inset of Fig. [IJa). Current experimentally con- 
firmed n-type donor and the direct STEM observation 
of Bis e i antisite defect supports a scenario that Bis e i 
might be viewed as a pair of Vsei and Bi interstitial, 
i.e., the Se vacancy acts as electron donor and Bi re- 
mains neutral. This scenario is reasonable, considering 
that the van der Waals radius of Bi (2.07A) is very close 
to the ionic radius of Se 2- (VI=1.98A).[15 Indeed, based 
on the formation energy calculated using first-principles 
calculations, Bis e i antisite defect has been shown to be 
favorable, compared to the Se vacancy defect with a Vsei 
character. [7 

The Bi2-layer intercalation must introduce severe lat- 
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tice distortion near regions where the van der Waals gap 
starts to dilate in order to accommodate the excess Bi, 
as shown in Fig. [2^c) . Bis e i antisite may introduce less 
strain to the system compared to the intercalated Bi in 
the van der Waals gap, but both can break the inver- 
sion symmetry locally, which might be associated with 
the observed symmetry-forbidden peaks found in Raman 
spectrum of thin film studies. [T6HT8] There are relatively 
more Bis e i aitisite defects found near the specimen edges 
than in the bulk from studies of multiple batches of either 
as-grown or 300° C annealed samples, as demonstrated in 
Fig. |5|b) . While surface defects could cause severe band 
bending away from the Dirac point at the surface, cur- 
rent findings based on STEM observation alone require 
other experimental techniques to verify. 

In summary, we conclude that crystals grown using 
Se-rich flux for the Bridgman growth method can reduce 
Se vacancy level significantly; however, it is not possible 
to prevent Se vacancy formation to achieve p-type dop- 
ing using excess Se under current the Bridgman growth 
condition. Two major types of defects are observed in 
Bi 2 Se3 crystals grown using Se-rich or Bi-rich flux: the 
randomly intercalated neutral Bi2-layer patches in the 
van der Waals gaps, and the Bi^i antisite point defects 
that act presumably as electron donor. Low temperature 
annealing of Se-rich flux grown crystal at 300° C helps to 
remove the intercalated Bi2-layer patches effectively and 
can create a large area of ordered near stoichiometric 
compound. 
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